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The nonlinear dispersion relation of geodesic acoustic modes 
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The energy input and frequency shift of geodesic acoustic modes (GAMs) due to turbulence in 
tokamak edge plasmas are investigated in numerical two-fluid turbulence studies. Surprisingly, the 
turbulent GAM dispersion relation is qualitatively equivalent to the linear GAM dispersion but can 
have drastically enhanced group velocities. In up-down asymmetric geometry the energy input due 
to turbulent transport may favor the excitation of GAMs with one particular sign of the radial phase 
velocity relative to the magnetic drifts and may lead to pulsed GAM activity. 
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Introduction. — Geodesic acoustic modes (GAMs), 
oscillating plasma flows associated with radial gradients 
of the electrostatic potential are an ubiquitous phe- 
nomenon in tokamak edge plasmas Q. By shearing the 
turbulence, GAMs are able to modulate the turbulent 
transport and to reduce the saturation intensity of the 
turbulence [3|. 

Using numerical two-fluid turbulence studies with the 
code NLET we investigate the energy input and fre- 
quency shift of GAMs due to their interaction with the 
turbulence. From the radial structure of global GAM 
eigenmodes Q in non-local (non-Boussinesq) computa- 
tions we deduce a nonlinear, i.e. turbulence induced, dis- 
persion relation. Although the observed dispersion agrees 
qualitatively with the ones obtained from linearized two- 
fluid equations [1, 01 , the corresponding group velocities 
tend to be surprisingly high. For realistic edge parame- 
ters, the global eigenmodes have a radial extent of several 
centimeters and could account for the frequency plateaus 
in ASDEX Upgrade Q , which deviate from the expected 
T^/^ dependence of acoustic frequencies. In turbulence 
runs with broken up-down symmetry we observe peri- 
odic bursts of turbulence and GAM activity. They can 
be traced back to an energy transfer from the turbulence 
to the flows which is odd in the radial wavenumber and 
absent with symmetry. This behavior resembles strik- 
ingly the experimentally observed pulsation during the 
I-phase in ASDEX Upgrade or the periodic turbu- 
lence suppression in NSTX [§]. 

Computational Setup. — We shortly summarize the 
basic computational setup of the turbulence runs. The 
ratio of the turbulence- to the background-gradient scale 
lengths is = L±/Ln oc p* , where L± is a characteris- 
tic turbulence scale length perpendicular to the magnetic 
field, = dr Inn, and the flux label r is the minor ra- 
dius at the outboard midplane. (A ratio of A ~ 1 — 50 is 
typical for the edge region.) The computational domain 
covers 490 sound gyro radii ps = {^mi{Te + T^))-^/^/ {eB) 
radially and poloidally. rm is the ion mass, Te,i the 
electron and ion temperatures, e the elementary charge, 
and B the magnetic field. The adiabatic exponent 7 



is 5/3. Furthermore, we choose r]i = Ln/Lxi = 2.4, 
rye EE Ln/Lr^ = 0, e„ = 2L„/i?o = 0.05, ad = 0.5, 
e-u = 0, and 7^ = 0, where ad is the ratio of the drift 
frequency to the ballooning growth rate, 7j, represents 
magnetic pumping, and £„ the sound speed (for defini- 
tions see (Jl). i?o is the major radius at the outboard 
midplane. The electrons are isothermal and treated non- 
adiabatically. Density and temperature are normalized 
to their values in the middle of the radial computa- 
tional domain, i.e. {n,f^/i) = X{n,T^/i)/{no,T^/i^o), and 
T = Ti_Q/Tefl = 1. The background density and tempera- 
ture profiles are linear, hh = A— 0.25r and Tn, = A— 0.62r. 
Thus, the local GAM frequency at zero radial wave num- 
ber ujGAM,oif) varies by roughly ±0.3wgaa/,o('' = 0), 
where r = marks the middle of the radial domain. The 
magnetic field is normalized to its value at the outboard 
midplane. With the time units to = Reff/{V2cs), where 
Cs = {"/{Te + Ti)/miY/'^ is the sound speed, l/i?e// = 
{2{C{9f)Y/'^/Ro,C{e) = (i?o|Vr|/S)(bxVlnB)-f and 
(...) = / • • • S~^d/|[/ § B~^dl\\ is the fiux-surface aver- 
age, ujgamA'^ = 0) = 1. 



Nonlinear GAM dispersion. — Figure [T] shows the re- 
sult of a nonlocal (A « 200) ITG turbulence computation 
with circular fiux-surfaces, which displays pronounced 
GAM activity. The radial variation of the GAM fre- 
quency caused by non-locality is obvious in the corre- 
sponding E X B-fiow profile [Fig. [Da)]. Its temporal 
Fourier transform [Fig. [T]b)] reveals that GAMs with fre- 
quency a;GAA/,o(''o) (indicated by the dashed line) radiate 
significantly outwards to r > rg but disappear for r < rp. 
The deviations from the local GAM frequency reach up to 
30 %. The application of a narrow band pass filter [Fig. 
[T] c)] shows that the fiow profile consists of global GAM 
eigenmodes. Such an eigenmode can be described by a 
WKB wave packet with a local wave number kr obeying 
i-^GAMir^kr) = const. To relate the observed deviations 
from ujGAMfl{To) to the nonlinear GAM properties, we 
choose the ansatz LOGAM{r, kr) = WGyiA/,o('^)(l -\-akf.) for 
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FIG. 1. (color online) Results of a non-local NLET turbulence 
studywith slow parameter variation, A « 200. a) Flux-surface 
averaged poloidal E x B velocity, b) Temporal Fourier trans- 
form (linear color scale) of a) with the local GAM frequency 
indicated by the dashed line, c) Band-pass filtered flow profile 
at oj = 0.98 with a fit to a curve of constant phase according 
to Eq. ^. 

the dispersion relation. The resulting solution for kr is 

where L^j is the gradient length of wgaa/.o- Depending 
on the sign of a, kr is imaginary on one side of the flux- 
surface r = ro and real on the other, kri^ro) = implies 
the existence of a reflection layer at tq . Using Eq. ([T]) and 
the condition f = uiGAM,o{ro)/kr ior the phase velocity 
yields the curve of constant phase 

ti^)-l\l-^ r . 2 , A r~r,f' + to. (2) 
ay aL^(ro)tJ^^M,o(^o) 

A least squares fit for the band-pass filtered flow profile 
in Fig. [T]c) with uJcAMfiifo) = 0.98 and the correspond- 
ing reference radius rg indicated by the dashed line yields 
a « 41p^, roughly a hundred times the value observed in 
the absence of turbulence (|a| p'l/2 Q). A comparison 
with the corresponding local run in Fig. [2] a) confirms 
the obtained dispersion relation. The above values for 
a can be used to estimate the radial extent of a global 



eigenmode whose radial structure is on one hand given 
by an Airy function Q. On the other hand the turbu- 
lence excites GAMs only in a limited range of radial wave 
numbers, whose distribution is assumed to be a Gaus- 
sian of width (Tfc centered at kr = ko- Since the Airy 
function Ai(x) drops only as x^^^, the wave number dis- 
tribution represents the dominant limitation of the mode 
width. The estimate resulting from Eq. ([T]) is therefore 
5r < AaLi^akko- On a basis of many turbulence com- 
putations we choose kops ~ 0.1 and akPs ~ 0.05. In 
tokamak edge plasmas L^i ^ Lt ^ 0.01 . . . O.li?. Thus, 
assuming a major radius of i? 1 m, we obtain mode 
widths of Sr < 0.1... 1mm, i.e. of the order of the 
gyro radius for a = plf^- In contrast, for a = 41 pi 
6r < 0.8 . . . 8.2 cm. The corresponding group velocities 
are of the order of the ion magnetic drifts and the dia- 
magnetic drift, respectively. Therefore, the width of non- 
linear eigenmodes can be comparable to the scale length 
of the frequency plateaus reported in [§] . 

Now, we develop an understanding of the turbulence 
induced frequency shift. The turbulent source terms are 
particularly transparent if we write the time evolution of 
the GAM in terms of a state vector 5* = {pgam,vgam) 
as 

dt^ = 5^ili,n + 5^„i, (3) 

which can be derived from the two-fluid system in Ref. 
[i]. Here, S^un = [vgam^-Pgam) describes the linear 
time evolution of the GAM and S'^ni = (sr + Sdia,su_) 
the nonlinear source terms. The GAM related part of 
the pressure fluctuations p is pgam = {C{9)p) with 
C = C/(C^)-^/^. The poloidal velocity component of the 
GAM is defined as vgari = {dr(p) = {vE,e/Cli where 
VE,e is the poloidal component of the E x B velocity 
= X V^(/). Here, is the potential fluctuation and 
^ = (1 + r)/(2(l + 7t)). Pressure and electric poten- 
tial perturbations, length, and time are normalized such 
that the GAM frequency logam.q = 1 and the free energy 
^ is E ~ p^ + [drcj))"^ ■ The divergence of the Reynolds 
stress is given by Sn = —{dr{vE,eVr,ion)), the divergence 
of the up-down asymmetric component of the turbulent 
transport by sr = ~{dr{G(9)vE,rP)), and the background 
diamagnetic velocity by Sdia = {ardripi + PrTi)) with 
Ur = r/(l + r) and Pt = 7r(l +7t) and the ion pressure 
Pi = 7i + Ti. The energy stored in the GAM corresponds 
to the squared length of the state vector VE". Obviously, 
S'^iin is orthogonal to 5" and is responsible for the GAM 
oscillation. Components of changing the length 

of 5* alter the GAM energy, the components orthogo- 
nal to ^E* cause frequency shifts. Using the basis vectors 
EE = (1/1*1)* and = {1/\^){vgam,-Pgam) one 
obtains the nonlinear energy input dtEni = (i5*„; • c-eY 
and the frequency shift Su = {{l/\^\)5'^ni ■e^)t with the 
time average (. . . )t- For the Fourier mode with w = 0.98 
[Fig. [T]c)], this splitting of the nonlinear terms into en- 
ergy inputs and frequency shifts clearly shows that the 



FIG. 2. (color online)a) Spectrum of the flux-surface aver- 
aged poloidal E X B velocity (logarithmic color scale) of a 
local NLET turbulence study otherwise equivalent to the one 
shown in Fig. [T] (using the parameters at r = 0) and disper- 
sion relations uJGAM{r, kr) = u!GAM,o{r){l-'rak'^) with a = pi 
(linear dispersion) and a = 'ilpi b) Nonlinear frequency shifts 
due to the couplings between turbulence and GAMs for the 
Fourier mode cu = 0.98 of Fig. [T]a): frequency shift necessary 
for the radial mode structure in Fig. [1] c) (dashed-dotted), 
shift due to the up-down asymmetric component of the tur- 
bulent transport (solid) , Reynolds stress by the (dotted) , and 
diamagnetic velocity (dashed). 



observed frequency shift is due to the coupling to the 
up-down asymmetric component of the turbulent trans- 
port sr alone, while the Reynolds stress and the diamag- 
netic velocity do not contribute [Fig. [2]b)]. The increase 
of the GAM frequency (and the phase velocity) due to 
the asymmetric transport F = VE,rP requires F to be in 
phase with pgam- A reduction of the resistivity results 
in reduced GAM group velocities, indicating a strong de- 
pendence of this effect on the ballooning character of the 
turbulence. 

Pulsed GAM activity. — We investigate the effect of 
up-down asymmetry of the magnetic geometry in local 
edge-turbulence computations using single-null geome- 
try with the same parameters as in the previous section 
but with adiabatic electrons. The flux surface is con- 
structed with the magnetic field of five toroidal current 
loops representing the plasma current and the currents in 
the shaping coils Q- The poloidal fiux is = ^ ■ = 
0.999*,ep, where A, = a,/2\n{{Ra + r,f + {Z + z,f) 
and Sep is the poloidal flux at the separatrix. The co- 
efficients Oi, Zi, and Vi measure the current of conduc- 
tor i, its position on the (vertical) Z-axis. and its ra- 
dial position. The radius of the magnetic axis is Ra- 
For consistency, the coefficients are chosen such that the 
forces on the plasma current exerted by the four shap- 
ing coils add to zero, i.e. J2i^i/i''~i + ^iY^^ = 0- We 
choose ro = zq ~ 0, gq = I, ri = 0, zi = to, ai = — m, 
r2 = -ra = Z2 = z^ = m/\/2, 02 = 03 = V2m, = 0, 
Z4 = —1, and 04 = 1 with m = 2 for lower (z^ —Zi 
for upper) X-point configuration. Ion curvature and V B 
drifts are directed in positive Z-direction. In local turbu- 
lence runs, the system saturates in both configurations, 
upper and lower X-point, by exciting GAMs. However, 
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FIG. 3. (color online) Flux-surface averaged poloidal E x B 
flows of NLET turbulence studies using the same parameters 
as in Fig. [T] but with a single-null geometry, a) Local run. 
Left: Ion magnetic inhomogeneity drift Vd directed away from 
the X-point. Right: Vd directed towards the X-point. b) 
Non-local run with A ~ 200 and towards the X-point. 
Color-coded: short-time- RMS of («(._§). Contours: turbulence 
intensity. Within the shaded areas, the turbulence intensity 
is above 70 % of its RMS average. 

a striking difference is the asymmetry of the nonlinear 
GAM drive with respect to the sign of the phase velocity 
of the GAM, which is obvious in the E x B-fiow pro- 
files in Fig. 13] a). Only GAMs with radially outward 
phase velocity are excited with the curvature drift di- 
rected away from the X-point whereas the phase velocity 
is radially inward with the curvature drift directed to- 
wards the X-point. The characteristic radial scale length 
of the GAMs is the same in both cases. The correspond- 
ing (fcr, a;)-spectra reveal that the GAM activity is con- 
centrated close to the linear dispersion. Since the linear 
group velocity of the GAM in single-null configuration is 
in direction of the ion curvature drift at a location oppo- 
site to the X-point , one can say that GAMs are excited 
such that the phase velocity is in direction of the group 
velocity. Including the parameter variation (A ~ 200) the 
GAM fiow amplitude and the turbulence intensity be- 
come correlated and pulsed at a frequency much smaller 
than the GAM frequency [Fig. [3]b)]. Detailed analysis 
of the poloidal flow shows that the characteristic radial 
scale fcr.i of the GAM in the initial phase {t < 100) is the 
same as in the local studies. The phase velocity is neg- 
ative. However, due to the frequency gradient the wave 
fronts of the flows are being tilted in the r — t-plane and 
kr decreases. In response, the GAM-turbulence equilib- 
rium readjusts itself towards transport levels even higher 
than in computations with artificially suppressed zonal 
flows. At fcr = the GAM and turbulence intensity 
reach a maximum and then drop for fc^ > as the 
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GAM is damped. After this burst, the described cy- 
cle repeats. The pulse frequency obviously depends on 
the GAM frequency gradient and on the initially growing 
wave number, i.e. cui, ~ 2Tr / {kr^iPs)drUJG am, oif)- A range 
of ujb ~ 0.01 . . .O.buJGAM seems to be readily realizable 
by adjusting drUJcAMfiif) and assuming kr^iPs ~ 0.1. 
For opposite sign of the ion curvature drift, the GAM 
wave number is restricted to kr > k^^i- In this case, the 
response of the turbulence intensity on the GAM wave 
number turns out to be weak, and the pulsing is subdom- 
inant. 

On the basis of Eq. ((3)), we investigate the depen- 
dence of the GAM-turbulence interaction on the sign of 
the phase velocity. Since the Reynolds stress turned out 
not to be important for the present parameters, we drop 
sn- Moreover, we approximate Sdia ~ ctrdrip) because 
the temperature term in the diamagnetic source term 
does not change the result qualitatively. The source 
term due to the turbulent transport F can be written 
as sr = ~dr{Cr). The transport itself is expressed as 
F = Fi F2, where Fi « fi(0)(l - Tdcdt)vaAM (for 
the Galilei invariance of this term see Ref. [3[), and 
F2 « f2(6')(l - Tdcdt)drVGAM [1, E^- Here, Tdc rep- 
resents the turbulence decorrelation time, which sets the 
time scale for the reaction of the turbulence to changes of 
the shear flow. The poloidal structure of F^ is contained 
in the function Ti{d). Furthermore, dt{p) ~ — 9r(F) and 
VGAM ~ ~iii^GAALo/^^GAM)PGAM, where ujgam,o is the 
linear, ujgam the (complex) nonlinear GAM frequency. 
Note that ^{logam) ~ ^^gama = 1 and "^{ujgam) > 
^(wG^Af) = Wi. Inserting the expressions given above 
into Eq. ([3]), the growth rate oipGAM evaluates to 

UJ, « i (Cf i) + klar (f 1 + TdcCf 2)) + 0(fc,^). 

(4) 

Only the term linear in fc^, which originates from sr, can 
be responsible for the asymmetry of the phase velocity 
of the turbulent GAM excitation. In up-down symmetric 
magnetic geometries C[6) is up-down antisymmetric (e.g. 
sin(6') / for circular flux-surfaces) . Since empirically Fi 
is approximately symmetric with respect to the outboard 
midplane of the tokamak and positive, (CFi) = and the 
related growth rate is zero. However, in single-null ge- 
ometry, C{6) becomes very small close to the X-point. 
Hence, (CFi) 7^ and an asymmetry in the phase veloc- 
ity can arise in Eq. ([4]). With b x V In _B (the direction of 
ion magnetic drifts) directed upwards - as is convention 
in NLET ~ and upper (lower) single- null geometry, G{9) is 
positive (negative) opposite to the X-point. Thus, (CFi) 
becomes positive (negative), and the sign of the asym- 
metric growth rate —{kr/2){CTi) agrees with the GAM 
properties observed in Fig. [3] a) . Note that the above 
calculation also yields a frequency shift due to the asym- 
metric transport Slog am = k'^{CT2) which is positive be- 
cause both, C{6) and r2(^), are negative (positive) above 
(below) the outboard midplane [10|. 



Summary and Conclusions. — Compared to linear 
predictions, numerical ITG turbulence studies display 
strongly enhanced GAM group velocities in case of non- 
adiabatic electron response with a dependence on the 
parallel resistivity. By analyzing the nonlinear frequency 
modification we could identify the mechanism responsi- 
ble for this nonlinear frequency shift. The interaction 
with the turbulence can raise the group velocity of the 
GAM from the order of the curvature drift velocity (^ 0.1 
km/s T - 100 eV, S ~ 1 T, i? ~ 1 m) in the hnear 
case 01 to the order of the diamagnetic drift velocity 
(^ 1 km/s), which is the typical scale of turbulent mo- 
tions. A global nonlinear eigenmode can be wide enough 
to form frequency plateaus as observed in ASDEX Up- 
grade [1] . The group velocity determines the flux of the 
energy stored in the GAM. Therefore, fast propagation 
might also play a role in considerations about the effi- 
ciency of externally driven GAMs as transport barriers. 

Up-down asymmetry of the magnetic configuration re- 
sults in an additional GAM growth rate causing a pref- 
erence for one particular sign of the phase velocity. In 
NLET runs using single-null geometry with the ion cur- 
vature drift directed upwards and upper (lower) X-point 
this interaction favors the excitation of GAMs with neg- 
ative (positive) group and phase velocities. Taking into 
account a radial frequency gradient, the phase velocity 
preference results in a pulsed activity of turbulence and 
GAMs for upper X-point. In contrast, for lower X-point 
pulsing is subdominant and the turbulence saturates into 
a quasi-steady state with the mean turbulence intensity 
being significantly lower than for upper X-point. The 
pulsed GAM activity reported here may be related to 
the pulsed GAM activity observed recently in ASDEX 
Upgrade Q. It might also be involved in the explana- 
tion of the quiet periods in NSTX Since the ion 
curvature drift was directed towards the X-point in all 
of the discharges analyzed in 0], it would be interesting 
to investigate the changes related to an inversion of the 
curvature drift in the experiment. 
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